Goh et al. Reply: We introduced in a recent Letter d the 
load distribution following a power law on scale-free (SF) net- 
works. In addition, it was conjectured that the load exponent 
8 is universal as long as the degree exponent 7 is in 2 < 7 < 3, 
based on real-world networks and in silico models. In the 
preceding Comment \2], Barthelemy argues that 8 is not uni- 
versal, sensitive to the details of SF networks. In this reply, 
we notice that the discrepancy is mainly caused by different 
usages of definition of load in 1 1] and |2]. Following the defi- 
nition used in |2], we agree with the result of [2], however, we 
find that the question of the universality of the load exponent 
is not settled yet. 

In 1 1 ], the load 1% of a vertex k includes N —\ packets leav- 
ing and another N — 1 packets arriving at the vertex, where 
N is the total number of vertices. However, those 2(N — 1) 
packets are not included in |2J]. While the difference of 
2(N — 1) can be neglected for vertices with large load i in the 
limit of N — » °°, however, in finite-size systems particularly 
those compatible with most real-world networks comprising 
N = 10 3 ~ 10 4 vertices, this difference could produce a differ- 
ent value of 8. We perform extensive numerical simulations 
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FIG. 1: The load exponent as a function of the mean number of edges 
(m) emanating from a new vertex for various degree exponents 7 in 
the BA model and different system sizes, N= 10 4 (COandA^ 10 5 
(■)• 


on a larger scale N = 5x 10 s than the size N =10 4 previously 
used in [1] for the static model with 7 w 2.5, following the def- 
inition in u||, and find that indeed 8 turns out to be lower than 
8 « 2.2 beyond error bar as argued in |2]. This behavior also 
occurs in the model introduced by Barabasi and Albert (BA) 
when the number of edges emanating from a newly added ver- 
tex is m > 2 in finite-size systems. However, we will show that 
the universal behavior of the load exponent is still likely as far 
as SF networks are sparse. 

The load exponent for SF tree has been obtained analyt- 
ically to be 8 — 2.0, independent of the degree exponent 

7 L3, yj. We investigate how the exponent value 8 = 2.0 
changes as the number of loops increases. We modify the 
BA model in such a way that a new vertex attaches one or 
two edges to existing network with probability 1 — p or p, re- 
spectively. The mean number of edges emanating from a new 
vertex is then (m) = 1+ p. We investigate how the load distri- 
bution changes as (m) varies. When p = 0, the network is tree, 
and the load exponent is confirmed to be 8 w 2.0. We find 
that 8 increases to 8 « 2.2 by increasing (m) to (m) w 1.1 
at which the edges connecting different branches of the tree 
structure form sparse loops in a nontrivial manner. The value 

8 « 2.2 turns out to be robust, independent of the degree ex- 
ponent 7 for 2 < 7 < 3. Such behavior persists as long as 
(m) is smaller than a 7-dependent critical value, (m) e , beyond 
which 8 depends on 7 as observed in |2|]. Moreover, we find 
that the plateau region of 8 w 2.2 is extended as the system 
size N increases as shown in Fig. 1. These data suggest that 
the universal behavior of 8 may hold in some finite region of 
parameter space in the thermodynamic limit, at least for the 
sparse BA model. Thus the possibility of the universal be- 
havior of the load exponent is still an open question. Further 
details will be published elsewhere |5]. 
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